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SUMMARY 

During certain phases of nuclear rocket operation, namely startup and 
shutdown, flow of hydrogen gas in the reactor core passages may change from 
turbulent to laminar flow as a result of fluid temperature increases in the 
passages. Experimental data required to validate assumptions used to estimate 
the pressure-drop characteristics of heated flow channels operating in (or 
near) the laminar -turbulent flow transition region appear to be lacking. 

This report presents results from an experimental investigation conducted 
to determine the pressure-drop characteristics of normal hydrogen gas flowing 
through an electrically heated test section operating in the transition region 
of laminar -turbulent flow. The experimental investigation was conducted at 
average operating heat flux values of 0.320, 0.485, and 0.640 Btu/( sec ) ( ft^ ) . 
The Reynolds numbers (based on inlet conditions) ranged from approximately 
600 to 5400, and the ratio of outlet to inlet fluid temperatures ranged from 
approximately 2.5 to 11.0. The test section used in these experiments was a 
Nichrome-V tube with an inside diameter of 0. 116 inch and a heated length of 
approximately 50 inches. 

From the experimental data, a steady- state curve of pressure drop as a 
function of inlet Reynolds number was obtained for each operating heat flux 
value. These curves are concave upward; the inlet Reynolds number correspond- 
ing to the point of minimum pressure drop ranged from approximately 1800 at 
the minimum operating heat flux of 0.320 Btu/( sec) (ft^) to 2300 at the maxi- 
mum heat flux of 0.640 Btu/( sec) (ft^ ) . 

A comparison was made between the experimental and analytically calculated 
values of pressure drop for several of the test runs. With the exception of a 
few points, the experimental and calculated pressure drops for all test runs 
agreed within ±10 percent . 



INTRODUCTION 


The efficient application of nuclear power for rocket propulsion requires 
a compact high-power- density reactor in which the propellant gas (hydrogen) is 
heated to extremely high temperatures. In the homogeneous reactor core proposed 
for the nuclear rocket, the heat generated in the solid fuel element array is 
removed "by the hydrogen propellant that flows axially through hundreds of small- 
diameter, parallel coolant passages in the core. A typical fuel element coolant 
passage design for the nuclear rocket core is depicted in reference 1. 

In a parallel flow system of the type described in reference 1, the static 
pressure drop across each flow channel in the core is set by the inlet and exit 
plenum pressures; that is, the static pressure drop across all flow channels 
must be identical and equal to the pressure difference from inlet to exit 
plenum. In such a system, nonuniform flow distributions among the hundreds of 
parallel coolant passages could result from a variety of conditions, (e.g., 
dimensional differences among the passages, nonuniform spatial heat generation 
rates, etc.) and could place a serious limit on the design performance of the 
nuclear propulsion system. 

During certain phases of nuclear rocket operation, flow of hydrogen gas 
through the parallel passages of the reactor core may lead to any one of sever- 
al distinct flow problems. For steady-state operation with no heat addition, 
the pressure drop in the passages is an increasing function of the flow rate 
in both the laminar and turbulent flow regimes; that is, the pressure-drop 
flow-rate curve has a positive slope in the laminar and turbulent flow regimes. 
For laminar flow with constant heat addition, however, the slope of the steady- 
state pressure drop - flow rate curve may be either positive or negative, 
depending on the value of the outlet-to-inlet fluid temperature ratio in the 
heated passage. The outlet-to-inlet fluid temperature ratio at which the slope 
of the laminar-flow pressure-drop flow-rate curve changes from positive to 
negative is herein termed the "critical fluid temperature ratio. " For values 
of outlet-to-inlet gas temperature ratio greater than critical, the steady- 
state pressure-drop flow-rate curve exhibits a negative resistance to flow 
(i.e., a negative slope). For values of the temperature ratio less than crit- 
ical, the slope of the steady-state pressure-drop flow-rate curve is positive. 
Hence, with a constant rate of heat input, the characteristic pressure-drop 
flow -rate curve with laminar flow is U-shaped (i.e., concave upward, with a 
minimum pressure drop occurring at the critical fluid temperature ratio) . 

Increases in the heat addition rate have been shown to result in other U- 
shaped pressure-drop flow-rate curves that exhibit the following characteris- 
tics: for a given pressure drop, the flow rate for the left leg of the U- 
shaped curve increases, and that for the right leg decreases as the heat flux 
is increased. If a small perturbation in heat flux is imposed, the flow rate 
for both legs tends to decrease. For the left leg, this tendency is in oppo- 
sition to the steady-state trend. These phenomena have been widely termed the 
"laminar flow instability problem" and will be discussed in more detail in the 
section GENERAL DISCUSSION OF PROBLEM. 
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Analytical studies of the laminar flow of gases through heated tubes have 
been made by numerous investigators, especially since the start of the devel- 
opment of nuclear rockets. Some of these studies are reported in references 2, 
3, and 4. In each case, the analytical results are in agreement with the dis- 
cussion presented in the preceding paragraph. An experimental study of the 
laminar flow of helium through a heated capillary tube (0. 012-in. i.d.) is 
reported in reference 5. With a constant electrical heat input to this capil- 
lary, the experimental data of reference 5 indicate the existence of two widely 
different laminar flow rates for a fixed value of pressure drop in the capil- 
lary. 


During nuclear rocket flight, two distinct operating periods must be ex- 
amined from the laminar instability viewpoint. The first, startup, may not be 
of significance if the period of low flow operation is traversed rapidly. The 
second, shutdown, may conceivably lead to serious difficulties, depending on 
the method of cooling employed to remove afterheat . One proposed scheme for 
removing the reactor afterheat involves pulse cooling; the use of pulses of 
appropriate strength and spacing may make it possible to maintain turbulent 
flow and, hence, stability. The cost of pulse cooling from a propellant stand- 
point depends directly on the mission intended for the nuclear rocket; whether 
or not this cost may be prohibitive cannot be stated at this time. 

A second method for removal of reactor afterheat may involve a continuous 
flow of coolant at an average rate which may be considerably smaller than that 
required for pulse cooling. If this scheme is employed, it is almost assured 

that at some time following shutdown mixed flow would exist in the reactor 

core; that is, as a result of fluid viscosity increases with increases in 
temperature, turbulent flow would enter the passage and laminar flow would 
leave the passage. 

A transition of this type (complete turbulent flow to complete laminar 
flow) would be expected to occur gradually over a length of flow passage 
equivalent to several tube diameters. Experimental data describing or defining 
the nature of a flow transition of this type appears to be lacking in the 
literature. Consequently, the number of tube diameters required for a complete 
flow transition to occur in a heated flow channel operating in the transition 

flow regime is uncertain. In addition, the value of the local bulk Reynolds 

number at which a flow transition of this type begins is also somewhat uncer- 
tain. 


In order to study the steady- state pressure-drop flow-rate characteristics 
of this type of flow in a passage the size of those currently used in the nu- 
clear rocket cores and for flow conditions anticipated during reactor after- 
heat removal, an analytical parametric study was conducted at the Lewis Re- 
search Center and is reported in reference 6. The study is dependent on an 
assumed laminar-turbulent flow transition criterion. For this study, flow was 
considered laminar at Reynolds numbers up to 2100 and turbulent down to about 
1000. In other words, a laminar flow remains laminar until the Reynolds 
number exceeds 2100, and turbulent flow remains turbulent until the friction 
factor in laminar flow is larger than that in turbulent flow. In the study 
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reported in reference 6, constant axial heat input was also assumed. The 
laminar-flow friction factor from the Poiseuille equation (ref. 7) and the 
Karman-Nikuradse formulation for use in turbulent flow with high film-to-bulk 
temperature ratios (ref. 7) were used in this analysis. Since the Reynolds 
number was based on film properties, wall temperatures had to be computed; 
these calculations were based on the assumption that the power generated in 
the solid is all transferred to the fluid. Furthermore, both laminar and 
turbulent heat -transfer correlations had to be employed. 

Calculations presented in reference 6 demonstrate that a minimum pressure 
drop always exists in a mixed flow passage. The value of the minimum pressure 
drop increases with heat input. 

In an effort to study in more detail the rather uncertain transition 
region between laminar and turbulent flow, that is, to determine the pressure- 
drop flow-rate characteristics and the relative position of the minimum steady- 
state pressure drop point, an experimental investigation was conducted at the 

Lewis Research Center with normal hydrogen gas flowing through an electrically 

heated 0.116 -inch-diameter Nichrome-V tube. This steady-state single-tube test 
marks only an initial step in a program of flow study through a nuclear rocket 

core. The purpose of this report is to present the experimental results of 

the initial phase of the overall program and compare these experimental results 
with theoretical predictions. The steady-state operating conditions covered 
in this first-phase investigation reported herein are as follows: inlet Reyn- 

olds numbers from 600 to 5400, outlet -to-inlet gas temperature ratios from 
2.5 to about 11.0, and average heat fluxes of 0.320, 0.485, and 0.640 
Btu/( sec) (ft 2 ). 

To determine reproducibility of the experimental results, tests for similar 
operating conditions were run three times. Results of the initial runs were 
reported in reference 8. Results from the two sets of reruns, together with 
the initial run results, are presented herein. 



Figure 1. - Static pressure drop in heated flow 
passage with laminar flow as function of flow 
rate for various relative heat-flux values. 


GENERAL DISCUSSION OF PROBLEM 
Steady State 

From an examination of the steady-state 
equations of static pres sure -drop for laminar 
flow in a uniformly heated flow passage, it can 
be shown that two different flow rates are pos- 
sible for a given value of pressure-drop in the 
flow passage. This phenomenon is illustrated in 
figure 1, where the pressure-drop in the flow 
passage AP is shown as a function of the cool- 
ant channel flow rate w for various relative 
heat flux values . (All symbols are defined in 
appendix A.) 

A derivation and general analysis of the 
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equation for static pressure drop of an ideal gas flowing in a constant area 
passage with constant, uniform axial heat input is presented in appendix B. 


The resulting steady-state expression for static pressure drop in a con- 
stant area passage with an assumed constant and uniform axial heat input 
(eq. (.1) of appendix B) is 
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The characteristic concave upward shape of the curves in figure 1 can be 
demonstrated from an analysis of equation 1. In appendix B, this equation is 
analyzed, and the following general conclusions are drawn: 


(l) With the heat input to the fluid approaching zero (i.e., with Q -* 0) 
such that the flow is essentially adiabatic, the static pressure drop AP S j.-2 
is an increasing function of flow rate w in both the laminar- and turbulent- 
flow regimes. 


(2) With a constant and uniform finite heat input to the fluid, however, 
the static pressure drop may be either an increasing or decreasing function of 
flow rate: 


(a) For the case in which the flow is turbulent throughout the pas- 
sage (i.e., for Reynolds numbers >2100) the static pressure drop AP _ _ 

• s , 1 — 2 

increases with flow rate w for all values of heat input. 

(b) For the case in which the flow is laminar throughout the passage 

(i.e., for Reynolds numbers <2100) the static pressure-drop AP _ 0 decreases 

s , 1 — 2 

with flow rate provided the heat input is sufficient to yield an outlet -to- 
inlet fluid temperature ratio which exceeds t 

crit 

The critical fluid temperature ratio ^ or laminar flow of normal 

hydrogen gas in a uniformly heated passage having a length and diameter approx- 
imately equal to a nuclear rocket core passage was determined by examining 
equation (l) for extreme values. A relative minimum pressure drop was 
found to correspond to a critical fluid temperature ratio of approximately 
4.0. The procedure used to calculate the critical fluid temperature ratio 
for laminar flow of normal hydrogen gas is described in appendix B. 

The analysis presented in appendix B is considered approximate in nature. 
As pointed out, the derivation of the static pressure-drop equation (eq. (l)) 
was based on several simplifying assumptions . It should also be recognized 
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that the computed critical fluid temperature ratio (of approx. 4.0) for all 
laminar flow was based on the assumption that the heat input was constant and 
uniform along the axis of the passage. It is expected that the value of t ^ 
would he somewhat different if the axial heat-flux distribution were nonuniform. 


Nonsteady State 

The transient facet of the laminar flow instability problem is that which 
applies to the stability of a heated flow channel to perturbations when oper- 
ating at a point on the steady-state pressure-drop flow-rate curve where the 
slope of the curve is negative. 

With reference to figure 1, it appears from the steady-state pressure- 
drop relations that if an operating point is chosen in the regime of negative 
slope (e.g., point A in fig. l) and the heat flux is increased slightly by an 
amount dQ, the flow rate w must increase if the pressure-drop AP is fixed. 
An analysis of the transient behavior of a heated flow channel does not., how- 
ever, lead to this conclusion. 

If AP is fixed (as would be the case in the multichannel parallel flow 
system with fixed pressure boundaries) and the heat flux in one channel of the 
parallel flow system is increased by an amount dQ, a decrease in w would be 
expected as a result of a decrease in the average fluid density in the flow 
passage . 

Results of unpublished analytical studies indicate that the negative slope 
region of the steady-state pressure-drop flow-rate curve is unstable. Positive 
perturbations in heat input have been shown to cause a continuous unbounded 
divergence in flow rate from the initial steady-state flow rate . At the low 
operating heat fluxes which are associated with the nuclear rocket afterheat 
removal period, the flow divergence caused by positive heat-flux perturbations 
occurs at a relatively slow rate and is accompanied by a continuously increas- 
ing flow passage surface temperature. The calculated unbounded flow-rate 
divergence associated with positive perturbations in heat input is based on 
the assumptions that the heat input to the passage is constant and that heat 
conduction from the passage is insignificant. 

The effect of heat interchange by conduction in the solid material of the 
reactor core would tend to prevent a continuing increase in temperature in 
those flow channels that are affected by an induced positive perturbation in 
heat input. 

The uncertainty in the conduction contact resistance between the solid 
fuel modules in the reactor core, however, makes it impossible to accurately 
determine the effectiveness of heat conduction in reducing the magnitude of 
nonuniformities in temperature and flow which might result from this phenomenon. 

From the previous discussion, it appears that laminar flow in the passages 
of the nuclear reactor core could present a serious problem, particularly during 
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the period of low power operation which follows a shutdown of the nuclear rock- 
et reactor. 

In the single -pas sage flow experiment described in this report., the 
laminar- flow perturbation instability problem was not encountered because of 
the relatively large flow impedance added into the system by the flow and pres- 
sure control valves. The pressure losses across these valves were relatively 
large compared with the pressure losses in the heated tube. As a result, the 
control valves had a stabilizing influence on the flow rate in the test system. 

Since the perturbation instability problem was not considered as an objec- 
tive in this experimental investigation, no further discussion of that facet 
of the problem will be made. 


EXPERIMEOTAL APPARATUS AM) PROCEDURE 
Test Assembly 

A sketch of the experimental setup is shown in figure 2. The hydrogen 
gas used in these experiments was supplied from a tube trailer, metered by 
calibrated rotameters, and precooled by a liquid nitrogen heat exchanger prior 
to entering the heated test section. The test section for these experiments 
was a Nichrome-V tube with an inside diameter of 0.116 inch, a wall thickness 
of 0.016 inch, and a heated length of 50 inches. The electrical resistivity 
of Nichrome-V at several temperatures from 70° to 2000° P is shown in table I. 
Prom an inspection of table I, it is seen that the resistivity of Nichrome-V 
is large and relatively constant over a wide range of temperatures. 

Two electrical power systems were used independently to supply power to 
the test section. In the initial set of test runs, the test section was heated 
by a manually controlled 6-kilovolt-ampere alternating-current power supply. 

For the two sets of reruns, a direct -current power supply with an automatic 
temperature feedback control was used. (it was mentioned previously that three 
separate sets of experimental data were taken for each operating heat flux. 

The primary purpose of running the two additional sets of data was to determine 
the reproducibility of the experimental measurements.) The direct -current 
power supply with automatic control was used to facilitate the setting up of 
the required power input for the reruns . The fact that two different types of 
electrical power were used (ac and dc) should have no effect on the measured 
data. 


The approximate positions of the electrical power connections are shown 
schematically in figure 2. To reduce heat sink effects, the electrical power 
tap on the downstream side of the heated test section was positioned just 
above the stainless-steel outlet plenum. Because of the large cross-sectional 
area of the outlet plenum, the electrical power dissipated in the plenum was 
relatively insignificant. (it was estimated that the power dissipated in the 
outlet plenum was less than 0.1 percent of the electrical power input.) The 
power tap at the test section inlet was located just above the liquid nitrogen 
heat exchanger and approximately 1.5 inches above the test section inlet. 
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Figure 2. - Schematic diagram of test assembly and instrumentation. 


TABLE I. - ELECTRICAL RESISTIVITY OF 
NICHROME-V AT SEVERAL TEMPERATURES 


[Data obtained from Driver -Harris Co. 
Technical Catalog, ECR58.] 


Temperature , 
°F 

Resistivity, 
( ohm) ( circular mi 

70 

650.0 

200 

660.5 

400 

674.0 

600 

685.0 

800 

693.0 

1000 

696.0 

1200 

692.0 

1400 

690.0 

1600 

693.0 

1800 

697.0 

2000 

701.0 


In order to minimize heat losses by free 
convection, the test assembly was enclosed in 
an evacuated chamber, which was maintained at 
an absolute pressure of approximately 10"^ 
millimeter of mercury. The heated fluid leav- 
ing the test section was cooled by a single- 
pass air-cooled heat exchanger and exhausted 
into a vacuum tank. The purpose of exhausting 
to the vacuum tank was to enable the testing 
to be conducted at a subatmospheric pressure 
level; this was considered necessary to max- 
imize the pressure drop across the heated test 
section. 

Figure 3 is a photograph of the test as- 
sembly showing the position of the Nichrome-V 
test section inside the vacuum enclosure. 
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Figure 3. - Photograph of test assembly. 


The relative locations and types of 
instrumentation used on the test assembly are 
indicated schematically in figure 2. The 
hydrogen flow rate was measured by two rotam- 
eters, each having a different range, con- 
nected in series and located in the flow line 
upstream of the test section. The rotameters 
were read remotely on a closed-circuit tele- 
vision receiver. 

The Nichrome-V test section was instru- 
mented with 36-gage (0 .005-in. -diam. ) iron- 
Constantan thermocouples, which were bonded 
to, and electrically insulated from, the out- 
side surface of the test section by a thin 
layer of ceramic cement. A total of 17 ther- 
mocouples, spaced 3 inches apart, were posi- 
tioned along the middle 48 inches of the 
50-inch-long heated length of tube. In order 
to minimize the conduction error in the wall 
temperature measurements, two turns of the 
insulated thermocouple lead wires were wrapped 
around and bonded to the surface of the heated 
tube at each measurement point. (During the 
testing, however, some of these thermocouples 
became detached or failed. As a result, wall 
temperature data were not obtained for each 
test run. ) 


The temperature of the hydrogen gas at the inlet and outlet plenums was 
measured by thermocouples positioned inside the plenums . A copper-Constantan 
thermocouple was used in the inlet plenum and an iron-Constantan thermocouple, 
surrounded with a stainless-steel radiation shield, was used in the outlet 
plenum. In addition, hydrogen gas temperature measurements were made with 
copper-Constantan thermocouples positioned in the transfer line immediately 
downstream of each rotameter. 


Static pressure taps were located in the plenums upstream and downstream 
of the test section. The pressure at the inlet plenum was measured with a 
transducer having a range of 12 pounds per square inch absolute (full-scale 
range); the pressure drop across the test section was measured with two differ- 
ent differential transducers having full-scale ranges of 0.7 and 0.3 pound per 
square inch. Pressure measurements were also taken in the transfer line imme- 
diately downstream of each rotameter. The electrical power input to the test 
section was determined from measurements of amperage and voltage drop across 
the power leads . 
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Precision of Experimental Measurements 


The recording instruments used in this experiment were of high precision. 
The pressure transducers ("both the differential and absolute transducers) were 
rated at ±1 percent of full scale. The accuracy of the thermocouples, partic- 
ularly those used to measure the test section wall temperatures, is difficult 
to determine. It was estimated, however, that errors caused by conduction 
through the thermocouple lead wires were relatively small. (As mentioned pre- 
viously, the thermocouple lead wires were wrapped around and bonded to the 
surface of the tube at each measurement point.) The rotameters used to measure 
the flow rates were calibrated with normal hydrogen gas; data points from the 
calibration were reproducible within about ±l/2 percent. It is estimated that 
the accuracy of the flow rate measurements were within ±3 percent . 

The reproducibility of the experimental measurement is not considered to 
be a quantitative indication of accuracy. Nevertheless, reproduction does 
provide a basis for confidence in the test results. For the most part, the 
reproducibility of the test data was reasonably good. 


Experimental Procedure 

For each steady-state test run, the static pressure and temperature at the 
test section inlet were maintained constant at approximately 10 pounds per 
square inch absolute and 140° R, respectively. The required flow rate for each 
run and the constant value of inlet pressure was established by opening the 
pressure and flow control valves to a proper setting. The electrical power in- 
put to the test section was then adjusted (either manually with the ac power 
supply or automatically with the dc power supply) to give an outlet fluid 
temperature for the established flow rate, which corresponded to a predetermined 
constant heat flux value. When the system reached steady state (i.e., when 
equilibrium temperatures and pressures were reached for the established flow 
rate) the experimental data were recorded. For the initial test runs, a con- 
tinuous 30-second analog recording of the data was made for each test condition. 
These data were taken on an automatic frequency -modulated tape recorder, proc- 
essed, and played back on strip charts by a pen-type oscillograph. The analog 
data acquisition and the play -back system are described in reference 9. For 
the two sets of reruns, the experimental pressure measurements were recorded on 
strip-chart recorders, and temperature measurements were displayed on a digital 
voltmeter . 

Experimental Investigations were made for three nominal values of operating 
heat flux (0.320, 0.485, and 0.640 Btu/(sec) (ft 2 ) ). For each test run, the in- 
let fluid temperature Tq was maintained constant (within ±1° R) at 140° R. 

The operating conditions for each of the three nominal heat flux values were 
determined prior to running the experiments. For each nominal heat flux, sev- 
eral values of flow rate were specified, and the corresponding values of outlet 
fluid temperatures (required to give the nominal values of heat flux) were de- 
termined from the following equation: 
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In running the experiments, some difficulty was encountered in establishing 
the exact values of the specified flow rates and corresponding outlet fluid 
temperatures. For some of the test runs ; the measured values of w and T^ 
were found to be slightly different from those computed prior to testing. 
Consequently, the actual Q/A s values (computed by use of equation (2) from 
the experimentally measured values of w and Tg) differed slightly from the 
reported nominal values. The actual Q/A s values for the individual test runs 
were determined and will be presented in the section RESULTS AM) DISCUSSION. 

Because of heat losses from the test section by radiation to the enclosure, 
and, to a lesser extent, by conduction from the inlet and outlet ends of the 
test section, the electrical power input could not be used as a measure of the 
average heat flux. 


PRESSURE -DROP ANALYSIS 

An estimate of the pressure drop In the heated passage was made for several 
of the experimental runs. Since the axial heat-flux distributions for the in- 
dividual test runs were apparently nonuniform, especially for those runs made 
at the higher operating temperatures, the fluid temperature profiles required 
for a pressure-drop analysis could not be computed directly from equation (2). 
(Fluid temperature measurements were made only at the inlet and outlet of the 
test section.) Therefore, the experimental wall temperature profiles (to be 
discussed in the section RESULTS AM) DISCUSSION) were used with appropriate 
convection heat-transfer relations to estimate the fluid temperature profiles 
and heat-flux profiles for the individual test runs. 

The estimated fluid temperature profiles were used in the analysis to 
compute the pressure drops across the heated test section for these test runs. 
The procedure used to calculate the pressure drops consisted of dividing the 
total test section length into 26 equal segments, each of length Ax, and com- 
puting the static pressure drop in each segment of the heated tube from inlet 
plenum to exit plenum. (The total length of the Nichrome-V tube, including 

the 1^-in. entrance length, was 51— in.; hence, Ax was approx. 2.0 in.) The 


pressure drops for the individual test runs were computed from the following 
equation: 
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The numerical operations required to compute the pressure drops from equa- 
tion (5) were performed on a high-speed digital computer. (Equation (3) is es- 
sentially the same as equation (B4) of appendix B. It is written in a somewhat 
different form herein as a convenience for use on the digital computing machine . ) 
The flow rates , inlet fluid temperatures and pressures used in equation (3) were 
obtained from the experimental measurements. 


For many of the test runs analyzed , the flow entering the heated tube was 
turbulent. As a result of fluid temperature increases along the tube length, 
the flow changed to laminar flow (based on local bulk Reynolds numbers) at some 
position downstream of the tube inlet . 


In the analysis, an assumed turbulent -to-laminar flow transition criterion 
was used. For the case of turbulent flow entering the test section, (i.e., for 


Re inlet 


>2100) the transition from turbulent to laminar flow in the heated pas- 
sage was assumed to occur at a position downstream of the passage inlet where 
the local bulk Reynolds number reached a value of 1000. For the case of laminar 
flow entering the test section^ i.e., for Re^ n -^ e ^ < 2100y the flow was consid- 
ered laminar throughout provided the local bulk Reynolds number at all points 
in the passage was less than or equal to 2100. These flow transition criteria 
are the same as those described and used in reference 6. As pointed out in 
reference 6, the local Reynolds number transition value of 1000 for turbulent - 
to -laminar flow was arbitrarily assumed. 


The friction factors used in the pressure -drop analysis were computed from 
the following relations: 


(ref. 


(l) For laminar flow with moderate wall-to-bulk temperature differences 

10 ) 


f = 


16 1 

\ ^wall 

Re, 

I , 


VI. 35 


(4) 


As stated in reference 10, the effects of heat addition together with 
hydraulic entrance region effects are combined in the temperature ratio term 
in equation (4). In appendix B, these combined effects on friction factor f 
were ignored, and the Poiseuille equation for isothermal laminar flow (eq. (B5)) 
was used. 


(2) For turbulent flow in long smooth tubes (ref. 7) 


f = 


0.0014 + 0,125 
„ 0.32 

RS b 


(5) 


Equation (5) is the Koo expression and is recommended for turbulent -flow 
Reynolds numbers ranging from 3X10 3 to 3X10 6 . Because of the assumed transition 
criteria, it was necessary (in some cases) to use equation (5) for Reynolds num- 
bers down to 10 5 . 

The Koo expression (eq. (5)) is somewhat different from the turbulent -flow 
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friction factor equation in appendix B. The Koo expression was used in this 
analysis primarily Because it more accurately describes the friction factor 
curve for smooth tubes at low turbulent Reynolds numbers. The hydraulic en- 
trance region effects with turbulent flow were determined from reference 11 to 
be relatively insignificant and were ignored in the analysis. 


For each of the test runs analyzed, the local heat-flux values and/or 
fluid temperature values required to compute the static pressure drops in each 
axial segment (of length Ax) of the heated passage were obtained from a simul- 
taneous solution of the following equations: 



*x 


m /^x+Ax/2 + ^x-Ax/ 2\ 

■hrall,x “ l g J 


Q_\ _ /_w_ 


(*./ 




■'"x+Ax/ 2 


C p dT 


"A -Ax/ 2 


( 6 ) 


(7) 


where x is the distance from the test section inlet . Values of the local 
heat-transfer coefficient h^ were computed from the following relations: 


(l) For laminar flow (ref. 12) 


hx ~ D- 


H 


0.036 


4 • 36 + 


(Re. Pr, i 
b,x b,x l 

x/D H 


1 + 0 . 001 ] 


7 Re Pr 

b,x b,x i 


x/D 


H 


( 8 ) 


Equation (8) was developed analytically and is based on an assumed constant 
axial heat input and a calculated nonparabolic fluid velocity profile. (As 
described in reference 13 the fully developed laminar flow velocity profile 
is not parabolic when heat is added. ) Equation (8) was used to compute 
local laminar-flow heat-transfer coefficients even though the assumption of 
a constant axial heat input (used in the development of eq. (8)) was not 
fulfilled in this experiment. 

(2) For turbulent flow with moderate wall-to-bulk temperature differences 
(ref. 7) 

h x = 0.023 ^ Re° -8 Pr°- 4 [l + ^°* 7 j (9) 

Equation (9) is recommended for complete turbulent flow (i.e., for Reyn- 
olds numbers > 10^; however, equation (9) was used in this analysis for Reyn- 
olds numbers as low as 10^. 
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Since values of heat-transfer coefficient in the transition flow regime 
have been found to vary considerably for different flow systems, the use of 
equation (9) in estimating turbulent heat-transfer coefficients at low turbu- 
lent Reynolds numbers is recognized as questionable. 

The assumed turbulent -to-laminar flow transition criterion used in the 
friction factor evaluation was also used in determining the appropriate equa- 
tion for computing the local heat-transfer coefficients. The values of wall 
temperature (^ va qq used in equation (6) to compute the local heat fluxes for 

the individual test runs were taken from the respective wall temperature curves. 

The local fluid temperatures computed from equation (7) were used in equa- 
tion (3) to determine the static pressure drop in each of the 26 axial flow 
segments; the computed static pressure drop from inlet to exit of the test 

section AP was obtained by summing the pressure drops in each of the 

s , 1 — 2 

flow segments. Pressure losses at the entrance and exit of the test section 
were included in the computed pressure drops and were estimated from the fol- 
lowing relations: 

(l) For the sudden contraction entrance pressure loss (from inlet plenum 
to tube) 


AP 


entrance 


0.5 


w 2 RT! 


144 2gP A 
s, 1 


(10) 


(2) For the sudden expansion exit pressure loss (from tube to exit plenum) 

( 11 ) 


, w 2 RTo 

AP . = 1.0 — 

exit 144 2gP A' 


'Z 


The pressure loss coefficients in equations (10) and (ll) were obtained from 
reference 14. 


Fluid properties used in the computations were taken from a digital- 
subroutine library of hydrogen properties. (The subroutine is described in 
ref. 15.) Local values of computed pressure and fluid temperature in the re- 
spective segments were used in the subroutine to evaluate the fluid properties 
used in the pressure drop and heat-transfer equations. 


RESULTS AND DISCUSSION 
Experimental Wall Temperature Data 

Figure 4 shows experimental wall temperature profiles for several test 
runs made at each of the three average operating heat fluxes. An inspection 
of the wall temperature profiles shows that the axial temperature distribution 
for each of the three heat fluxes becomes increasingly nonlinear at the higher 
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Temperature, 




(b) Average heat flux, ~0. 485 Btu/(sec)(ftA 
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Fractional distance from upstream power tap, x/L 

(c) Average heat flux, ~0. 640 Btu/($ec)(ft 2 ). 


Figure 4. - Test-section wall temperature profiles for different 
values of inlet Reynolds number. 


operating temperatures. Figure 4 shows 
that at the higher operating temperatures 
the slopes of the wall temperature curves 
are steep near the test-section inlet and 
decrease along the length of the heated 
tube . 

Since the local slopes of the wall 
temperature curves are a relative indica- 
tion of the local axial heat flux, it 
appears that for those test runs made at 
high operating temperatures (and/or low 
inlet Reynolds numbers) a major portion of 
the total heat input to the fluid was 
added in the leading half of the test sec- 
tion. The increasing nonlinearity of the 
wall temperature profiles at the higher 
operating temperatures can be attributed 
primarily to appreciable heat losses by 
radiation from the downstream portion of 
the tube surface. Therefore, even though 
the integrated heat input (as defined by 
eq_. (2)) was essentially constant for each 
test run in a particular heat-flux series, 
the distribution of heat input to the fluid 
would be expected to be somewhat different . 


Axial Heat -Flux Distributions 

Figure 5 shows the calculated axial 
heat -flux distributions for three of the 
test runs made at an average operating 
heat flux of 0.640 Btu/( sec) (ft^) . Al- 
though the areas under each of the heat- 
flux curves in figure 5 are approximately 
the same, the distributions are seen to be 
vastly different. The increasing nonuni- 
formity in heat flux with increasing oper- 
ating temperatures (or decreasing inlet 
Reynolds numbers) is attributed primarily 
to significant radiation heat losses from 
the downstream end of the heated test sec- 
tion. 

An examination of calculated heat- 
flux distributions for the other average 
operating heat fluxes indicates that the 
distributions in figure 5 are typical; 
that is, the calculated distributions of 
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Figure 5. - Calculated heat-flux distributions along length of test section for different values of in- 
let Reynolds number. Average heat flux, ~0. 640 Btu/tsecXft^). 


heat input for each operating heat flux become increasingly nonuniform with in- 
creases in operating temperature. 

Experimental Pressure-Drop Flow-Rate Characteristics 

The results of the laminar-flow instability experiments are presented in 
figure 6. The upper curves in figure 6 show the static pressure drop across 
the heated test section (expressed as the product of inlet pressure and static 
pressure drop) as a function of inlet Reynolds number for average heat-flux 
values of 0.320, 0.485, and 0.640 Btu/(sec) (ft 2 ) , respectively. An inspection 
of the pressure-drop curves shows that the position of the relative minimum 
pressure drop moves upward and to the right with increasing heat flux. At the 
minimum operating heat flux, the relative minimum pressure drop corresponds to 
an inlet Reynolds number of approximately 1800; at the maximum operating heat 
flux, the relative minimum pressure drop corresponds to an inlet Reynolds num- 
ber of approximately 2300. 

The lower curves in figure 6 show the experimental outlet -to -inlet fluid 
temperature ratios as a function of inlet Reynolds number. An examination of 
these curves indicates that the fluid temperature ratio Tg/T^ at the point of 
of minimum pressure drop increases from a value of 4.0 at the minimum heat flux 
to approximately 5.5 at the maximum heat flux. 

From the analysis presented in appendix B, the value of Tg/T^ correspond- 
ing to the point of minimum pressure drop for all laminar flow was estimated to 
be approximately 4.0. 
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TABLE II. - TEST DATA FROM HYDROGEN GAS PLOW INSTABILITY EXPERIMENTS 


Ru 

n Measured 

Calculated 

Measured 


flow rate 

average 

inlet 



heat flux. 

pressure. 


Ib/sec 

Q/A s , 

P S,1' 



Btu/(sec)(ft 2 

) psia 

1 

10.56x10-* 

0.330 

10.03 

2 

13.45 

. 320 

10.00 

3 

15.95 

.314 

9.86 

4 

19.70 

.330 

9.87 

5 

23.90 

.311 

9. 93 

6 

27.90 

. 315 

10.00 

7 

33.00 

.319 

9.90 

a 

41.05 

.311 

9.90 

9 

47.80 

.315 

10.00 

10 

53.50 

.322 

9.94 

11 

61.00 

.329 

10.00 

12 

10.21 

.327 

9.87 

13 

12.98 

. 329 

9-97 

14 

15.90 

.319 

9.90 

15 

19.84 

.317 

10.17 

16 

23.83 

.319 

10.12 

17 

27.80 

.318 

10.05 

18 

21.80 

.320 

9.87 

19 

17.50 

. 320 

10.05 

20 

11.35 

.320 

9.90 

21 

33.80 

.328 

10.05 

22 

64.70 

.311 

10.00 

23 

51.70 

.312 

10.00 

24 

46.60 

.314 

10.00 

25 

40.20 

.311 

9.99 

26 

33.90 

.315 

10.03 

27 

13.86 

.488 

9.96 

28 

14.96 

.477 

9.87 

29 

16.17 

. 471 

9.87 

30 

20.91 

.472 

9.87 

31 

25.80 

.487 

9. 96 

32 

29.82 

.481 

9.99 

33 

34.35 

.482 

9.96 

34 

38.50 

.480 

9.88 

35 

43.40 

.489 

9.96 

36 

51.80 

.484 

10.00 

37 

59.80 

.479 

9.85 

38 

75.80 

.488 

9. 97 

39 

17.61 

. 491 

10.00 

40 

20.85 

.480 

10. 19 

41 

25.82 

.483 

9.99 

42 

29.81 

. 482 

9.90 

43 

34.42 

.485 

9.90 

44 

36. /2 

.490 

10.00 

45 

46.25 

.487 

10.00 

46 

55.22 

. 482 

10.00 

47 

6 7.70 

.473 

10.00 

48 

82.00 

.484 

10.00 

49 

13.35 

.481 

10.03 

50 

15.40 

.483 

9.98 

51 

16.00 

.485 

9.90 

52 

16.86 

. 627 

10.10 

53 

19.20 

.640 

10.00 

54 

22.15 

. 638 

10.00 

55 

26.97 

.626 

9.94 

56 

32.80 

.630 

9. 92 

57 

38.80 

.626 

10.00 

58 

45.20 

.636 

10.00 

59 

93.40 

.653 

9. 96 

60 

21.60 

.646 

10.00 

61 

25.60 

.651 

10.00 

62 

31.10 

.644 

10.00 

63 

39.20 

. 646 

9. 98 

64 

43.70 

.647 

10.00 

65 . 

54.20 

.643 

10.00 

66 i 

55.20 

.639 

10.00 

67 

76.00 

. 642 

9. 96 

68 

72.20 

.642 

10.00 

69 1 

32.00 

.659 

10.00 

70 : 

37.50 

.654 

9.82 

7i : 

16.80 

.631 

9.90 

72 : 

18.20 

.633 

9.93 

73 < 

30.00 

.642 

9.96 


Measure^ 

d Calculate' 

d Measure' 

d Calculate< 

i Calculated 

static 

inlet 

outlet 

fluid 

outlet 

pressur 

e Reynolds 

fluid 

temper- 

Reynolds 

drop, 

AP s,1-2 

psi 

number 

temper- 

ature, 

t 2 , 

°R 

ature 

ratio, 

t 2 / Ti 

number 

0.410 

605 

1279 

9.15 

129 

.258 

7 70 

1036 

7.40 

188 

.220 

915 

900 

6.43 

245 

.200 

1130 

7 94 

5.68 

290 

.173 

1370 

657 

4.70 

457 

.166 

1600 

589 

4.21 

575 

. 162 

1890 

534 

3.82 

725 

.170 

2400 

449 

3.21 

1022 

.182 

2740 

420 

3.00 

1235 

. 195 

3070 

399 

2.85 

1441 

.206 

3500 

361 

2.58 

1750 

.350 

586 

| 1351 

9.65 

122 

. 244 

742 

I 1077 

7.67 

170 

.199 

911 

913 

6.52 

242 

.171 

1137 

767 

5.48 

345 

.163 

1365 

670 

4.78 

449 

.158 

1592 

600 

4.28 

564 

.175 

1250 

718 

5.12 

389 

.196 

1003 

850 

6.08 

339 

.305 

650 

1211 

8.66 

138 

.160 

1937 

535 

3.82 

730 

.232 

3710 

349 

2.49 

1900 

.191 

2960 

397 

2.83 

1390 

. 174 

2670 

424 

3.03 

1190 

.164 

2304 

463 

3.31 

972 

.152 

1943 

519 

3.71 

757 

.445 

7 95 

1460 

10.43 

157 

.401 

858 

1318 

9.41 

179 

.363 

927 

1243 

8.88 

201 

.289 

1200 

1004 

7.17 

300 

.260 

1480 

869 

6.21 

406 

.231 

1710 

770 

5.50 

508 

.222 

1970 

694 

4.96 

624 

.229 

2210 

636 

4.54 

754 

.237 

2490 

591 

4.22 

694 

.252 

2970 

521 

3.72 

1158 

.274 

3430 

472 

3.37 

1420 

.311 

4350 

412 

2.94 

2001 

.36 7 

1009 

1198 

8.56 

231 

.273 

1195 

1023 

7.31 

292 

.245 

1480 

865 

6.18 

406 

.230 

1708 

772 

5.52 

510 

. 238 

1973 

696 

4.97 

631 

.226 

2104 

669 

4.78 

690 

.251 

2650 

566 

4.04 

983 

.276 

3165 

498 

3.56 

1276 

.314 

3880 

434 

3.09 

1716 

.360 

4 700 

392 

2. 80 

2210 

.532 

764 

1495 

10.68 

150 

.420 

882 

1325 

9.47 

184 

.371 

917 

1289 

9.21 

195 

.564 

967 

1530 

10. 93 

188 

.465 

1100 

1384 

9.89 

225 

. 390 

1270 

1204 

8.60 

281 

.347 

1546 

1023 

7.31 

382 

. 310 

1880 

882 

6.30 

504 

. 295 

2225 

768 

5.49 

664 

.300 

2590 

695 

4.97 

827 

.481 

5360 

433 

3.09 

2365 

.445 

1238 

1273 

9.09 

263 

.357 

1467 

1106 

7.90 

334 

. 329 

1782 

938 

6. 70 

455 

.319 

2247 

785 

5.61 

662 

.318 

2505 

722 

5.16 

766 

.338 

3105 

615 

4.39 

1083 

.378 

3740 

539 

3.82 

1410 

.430 

4355 

488 

3.49 

1758 

.435 

4140 

503 

3.40 

1705 

. 370 

3550 

565 

4.04 

1295 

.335 

2150 

821 

5. 87 

618 

.657 

962 

1552 

11.08 

1B4 

.548 

1042 

1450 

10.35 I 

208 

.535 

5160 

439 

3.13 

2260 


Measured inlet 
pressure times 
measured static 
pressure drop, 
P s,l(A p s,l- 2 ) 


4.12 

2. 58 
2.17 

1.97 

1.72 

1.66 

1.60 

1.68 

1.82 

1.94 

2.06 

3.46 

2.43 

1.97 
1.74 

1.65 

1.59 

1.73 

1.97 
3.02 

1.61 

2.32 

1.91 

1.74 
1.63 

1.53 

4.43 
3.96 

3.58 
2.85 

2.59 
2.31 
2.21 
2.26 
2.36 

2.52 

2.70 
3.10 

3.67 

2.78 

2.45 
2.28 
2.36 
2.26 
2.51 

2.76 

3.14 

3.60 

5.34 
4.19 

3.68 

5.70 
4.65 
3.90 

3.45 

3.08 

2.95 

3.00 

4.79 

4.45 

3.57 

3.29 

3.18 

3.18 

3.38 

3.78 

4.28 

4.35 

3.70 

3.29 

6.50 

5.44 

5.33 
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Figure 6. - Experimental static pressure drop and experimental fluid temper- 
ature ratio against inlet Reynolds number with heat flux as a parameter. 


In figure 6, no distinction is made between the data points obtained from 
the initial set of test runs and those obtained from the two sets of reruns. 

In general, no significant difference was observed between the initial test run 
data and that obtained from the reruns; hence, it appears that the test points 
are reasonably reproducible . 

A summary of the measured and calculated data for each test run is presented 
in table II. An inspection of these data shows the experimental heat-flux value 
for each test run to be within ±3 percent of the reported nominal values. The 
procedure used in calculating the values in table II is demonstrated in appen- 
dix C by a sample calculation. 
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Figure 7. - Comparison of experimental and cal- 
culated pressure-drop values for three oper- 
ating heat fluxes. 


Figure 8. - Comparison of experimental and cal- 
culated fluid temperatures at test-section exit 
for three operating heat fluxes. 


Comparison of Experimental and Calculated Pressure Drops 

The calculated pressure drops for several test runs made at each of the 
three operating heat fluxes were compared with the respective measured pressure- 
drop values; the results of this comparison are shown in figure 7. The assump- 
tions^ equations, and analytical procedures employed in the calculations are 
described in the section PRESSURE-DROP ANALYSIS. As mentioned in the section 
EXPERIMENTAL APPARATUS AND PROCEDURE, wall temperature measurements were not 
obtained for all of the experimental runs because of thermocouple failures. As 
a result, pressure drops for all test runs could not be computed. 

Figure 7 shows that with the exception of a few points, the calculated 
pressure-drop values agree with the experimental values within about ±10 percent. 
The differences between the experimental and calculated pressure drops for the 
three different operating heat fluxes appear to follow no consistent pattern. 

In these experiments, no measurements were taken of local fluid temper- 
atures along the length of the heated passage. As a result, the calculated 
fluid temperature profiles and heat-flux profiles could not be experimentally 
validated. Fluid temperature measurements taken at the exit of the test section 
did, however, compare rather closely with the calculated outlet fluid tempera- 
tures . 

Figure 8 shows the differences between the calculated and measured values 
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of exit fluid temperatures for the test runs analyzed. From figure 8, it is 
seen that with the exception of one of the test runs ; the experimental and cal- 
culated fluid temperatures at the test section exit agree within ±15° R. 

In judging the comparison between experimental and calculated pressure 
drops in figure 7 , it should be recalled that the computed pressure-drop values 
were dependent on the heat -transfer relations used as well as the pressure drop 
relations. 

Nevertheless; the experimental and calculated values of pressure drop for 
the individual test runs are reasonably close. The data in figure 7 indicate 
that the pressure-drop and heat-transfer relations along with the assumed flow 
transition criterion used in the analysis were reasonable. 


DISCUSSION OF TEST DATA 

Some idea of the large flow and temperature nonuniformities possible in 
a system of parallel flow passages can be seen from an examination of figure 6. 
For example; the curve constructed through the data points for the average heat 
flux of 0.640 Btu/(sec) (ft^) shows that for a value of P AP 9 of 4.5; the 

flow rates and/or inlet Reynolds numbers that satisfy these conditions may 
differ by a factor of approximately four. The corresponding outlet fluid tem- 
peratures for the two flow rates differ by approximately 900° R. 

The large flow and temperature differences cited in the preceding example 
are based on the assumption that each of the parallel flow passages is thermal- 
ly isolated from the other flow passages in the system. In the solid core of 
the nuclear rocket; the magnitude of the nonuniformities in flow and temper- 
ature would be reduced by heat conduction between adjacent fuel cells. If the 
heat interchange by conduction between fuel cells is large enough to prevent 
thermal damage in the core; the problem of laminar-flow instability (as applied 
to the solid core of the nuclear rocket) may be of little significance. The 
uncertainty in the conduction contact resistance between the fuel cells in the 
reactor core makes it impossible to estimate accurately the effectiveness of 
radial heat conduction in preventing temperature and flow maldistributions of 
the type associated with the laminar instability phenomenon. 

Therefore; the magnitude of possible nonuniformities cited in the preced- 
ing example should be carefully interpreted. Application of these data and 
conclusions drawn regarding possible nonunif ormities of temperature and flow 
among parallel flow passages is justified only where the effect of thermal 
coupling between channels is negligible. 


CONCLUDING REMARKS 

Since gas-cooled nuclear reactors are normally designed for operation in 
the turbulent flow regime; the laminar -flow instability phenomenon reported 
herein would not appear to present a problem for steady-state design operation. 
During the startup and shutdown of a gas -cooled reactor; the system flow rate 
must pass through the laminar flow regime. In the case of shutdown; the re- 
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actor decay heat may be sufficiently large to require cooling of the core as- 
sembly for relatively long periods of time. If coolant consumption is to be 
minimized during the period of afterheat removal (as would be required for the 
missions being considered for the nuclear rocket) , it is necessary that the 
flow rate be controlled to a minimum value which would maintain the core assem- 
bly slightly below the maximum allowable temperature. 

The experimental data presented in this report show that a heated flow 
channel operating at steady state with a constant heat flux may have two widely 
different flow rates which yield the same value of pressure drop. As a result, 
it seems reasonable to expect that nonuniformities of temperature and flow 
could occur in a system of uniformly heated parallel flow passages operating at 
low flow rates between fixed pressure boundaries. 

For the range of test conditions investigated, the experimental data indi- 
cate that 

1. The slopes of the curves of pressure drop against flow rate and/or 
inlet Reynolds number are negative for outlet -to -inlet fluid temperature ratios 
greater than 5.5. 

2. The flow rate (or inlet Reynolds number) corresponding to the point of 
minimum pressure drop increases with increasing values of heat flux. At the 
lowest heat flux, the relative minimum pressure drop occurred at an inlet 
Reynolds number of approximately 1800; at the highest heat flux, the relative 
minimum pressure drop occurred at an inlet Reynolds number of approximately 
2300 . 


Several test runs were made at each average operating heat flux to inves- 
tigate the reproducibility of the experimental measurements. In general, the 
experimental pressure -drop measurements were reproducible, even in the regime 
of mixed flow. 

A comparison was made of experimental and calculated values of pressure 
drop in the heated flow passage. Fluid temperature profiles required for the 
pressure -drop analysis were estimated from a convection heat -transfer analysis 
in which measured values of the wall temperature distributions were used. An 
assumed turbulent -to-laminar flow transition criterion was used in the pressure- 
drop heat-transfer analysis to compute axial fluid temperature distributions 
for those test runs having mixed flow. The agreement in experimental and 
calculated values of pressure drop indicates that the heat-transfer and 
pressure-drop relations along with the assumed flow transition criterion were 
reasonable. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, November 23, 1965. 
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APPENDIX A 
SYMBOLS 

flow area ; ft 2 

inside surface area of test section, ft 2 

specific heat of gas at constant pressure, Btu/(lb)(°R) 

hydraulic diameter of test section, ft 

friction factor 

constant in friction factor equation (eq. (B5)) 
gravitational constant, ft/ (sec) 2 
enthalpy, Btu/lb 

local heat-transfer coefficient, Btu/(sec) (ft 2 ) (°R) 

thermal conductivity of gas, Btu/( sec) (ft ) ( °R) 

heat -transfer length of test section, ft 

exponent on temperature in viscosity relation (eq* (B6)) 

exponent on Reynolds number in friction factor relation (eq. (B5 ) ) 

static pressure, psia (except in appendix B where psfa is used) 

Prandtl number, pC^/h 

rate of heat transfer to gas, Btu/sec 

specific gas law constant, ft/°R 

Reynolds number, wD^/Ap 

temperature, °R 

fluid velocity, ft/sec 

weight flow rate, lb/sec 

axial distance from inlet of test section, ft 
viscosity, lb/ (ft) ( sec) 


\i Q constant coefficient in viscosity relation (eg.. (B6)) 
t outlet -to-inlet fluid temperature ratio, 

Subscripts : 
ave averaige 

Id Bulk (when applied to fluid properties, indicates evaluation based on 

bulk temperature) 

crit critical 

f film (when applied to fluid properties, indicates evaluation based on 

film temperature , T f = ( T wa ll * T t))/ 2 ) 

wall inside surface of test section 

x local value at distance x from test -section inlet 

1 test-section inlet 

2 test-section outlet 
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APPENDIX B 


APPROXIMATE CLOSED -FORM PRESSURE DROP -FLOW ANALYSIS 


Consider a heated flow channel of length L in which a fluid flowing at 
a rate w enters the channel at temperature Tp and pressure Pj_ and 
leaves the channel at temperature Tg and pressure Pg as shown in figure 9. 



If changes in potential energy of the fluid 
are neglected, the static pressure drop across 
the heated flow channel AP is given 

S , _L — c* 

by the following equation: 

AP = AP ^ + AP . _ (Bl) 

s,l-2 momentum friction 

If the fluid is assumed to behave as a per- 
fect gas and if the pressure -drop is small 
relative to the inlet pressure. 


AP 


momentum 


/.\2 * 

2 m £ i (To - 



(B2) 


AP 


friction D \A/ g P 
Hence, equation (Bl) may "be written as 


/.\2 /»x=L 

4 /w\ R 1 / ^ 

W gP, + / 

1 ^ *x=0 


fT dx 


AP 

s 


> 


1-2 





(B3) 


(B4) 


The friction factor f in equations (B3) and (B4) is expressed as a 
function of the local Reynolds number Re v by 


f = 



For laminar flow (ref. l), 



f = constant = 16 
n = 1.0 


(B5) 
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For turbulent flow (ref. 7) 


f Q = 0.046 
n = 0.2 


The local fluid viscosity p x is represented as a function of the local fluid 
temperature T by 

(® 6 ) 

Combining equations (B4), (B5), and (B6) yields 


AP 


Sjl-2 



R 1 

« P ! + P 2 


T„ 


T 


+ 


2 




(B7) 


If the Maeh number in the heated channel is small relative to one, the static 
temperature rise of the fluid is related to the heat input hy 



(B8) 


Integrating equation (B8) (assuming CL constant) and rearranging give 


m m 

2 1 wC 


(B9) 


P 


(For normal hydrogen gas, the assumption that C p is constant is rather poor 
at temperatures belov about 500° R. ) 


For a constant diameter channel, the surface area for heat transfer Is 


A = JtDL 
s H 


(B10) 


Multiplying numerator and denominator of the right 
A s = jcDj^L gives 

QtcDttL 

T - T = n 
2' 1 A ¥C 

s p 


side of equation (B9) by 


(Bll) 


The flow area of a circular channel is A = (it/ 4 ) (DjjLj hence, 

it = 4A/(D h ) 2 

Substitution of 4A/(D ) 2 for it in equation (Bll) gives 

II 
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(B12a) 


T 2 - T 1 " 


4QL 


A C Dj- 
s p HU 


If the heat flux is defined as Q/A s and the heat input is assumed to 
uniformly distributed along the length of the channel (such that 

(Q/A s ) = (Q/A^) \ equation (B12a) may be written as'' 
ave x/ 


T - Tp = 


4Qx 


v p m! 


The integral in equation (B7) is evaluated as follows: 

// x=L 

x=0 


T 1 + 


4Qx 




Mn+1 


dx 


/* X L ^Mn+l 
x=0 


dx = 


A s c p% A /rn J. 4Qx 


vMn+2 


Mn + 2 4Q, 


I T ± + 


ACL -1 
s p H a/ 


x=L 


x=0 


/ x=L 
0 


1 A s C p D H I 


T Mn+1 dx = 0 

x=0 Mn + 2 4Q 


, T + 4QL \ Mn+2 Mn+2 

lT l A C Dj*\) ' Tl 

S P h^a )J 


Combining equations (B12a) and (B13c) yields 


/ X=L T^ 1 d* = - ,(lf” +2 - if +2 ) 

y x=0 (Mn + 2)(T 2 - T 1 ) V ^ 1 / 

Substituting the result of equation (B13d) into equation (B7) gives 

2 


AP — 2 — 1 

^8,1-2 " * 


w 


S P 1 + P 2 ' A 


T 2 - T 1 


2 o^o 


d h(-d h ' (Mn + 2)(T ? - T i )^ 2 

L A 


rpMn+2 _ rpMn+2 


be 


(BIZb) 


(B13a) 


(B13b) 


(B13c) 


(B13d) 


(B14) 


26 



Multiplying numerator and denominator of equation (B14) by 


AP 


= pR T 1 /vfi T 2- T l 

h l-2 g P + P \ A I ) T 

-L C t J- 






:(! D H ) n ^ Mn + 2 ) t i ( t 2 - VlS T i > 


Tc 


iMn+2 


^ - 1 


From equation (B6), p = p T^; hence, 

X o 


1 o 1 


(B15a) 


(B15h) 


Substituting the result of (B15b) into (B15a) gives 


R T 1 


/ -\2 \T 9 - T, 

AP q n o = 2 ±i ± ( 2) )-£ i 

s ,l-2 g P + P \Aj \ T 


’1 2 


2f o^ LT l 


D h(I + 2)(T 2 - T x ) 


n 


/m \ Mn+ E 

- 1 

V T i/ 


Combining equation (B12a) with (B16) 


AP 


2R T 1 




w 


4QL 


2 f o |J 'l L A s C p D H T l a 


s,l-2 g p i + Pg \Ay ) A C D T E D h 

( S p HI A 


(! °n) n 


(Mn + 2) 


4QL 


(B16) 


X 


1 + 


4QL 


vMn+2 


A C D T - 

s p H 1 a 


- 1 


(1) 




Equation (l) is the general equation of the family of curves shown in figure 1. 


Taking the partial derivative of AP s 1-2 with respect to w/A at a 
constant heat flux shows that a relative minimum exists for the curves in fig- 
ure 1. 
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a ap 


s jl-2 


Hi 


= 2 5 T 1 


4QL 


Q/A c 


S P x + P 2 1 A„D tt C_T 


2f o^l M s D H C p T l^) 

8 D *) n 


XI 1 + 


s H"p~l 4 QLD h 

4QL • Mn+1| “ 


A G D T 1 ~ 

spVl A / L "sVh^iVa 

■ .\ 2”| 


4QL 


A C D T I t'' 2 


, Wig 

4QLD (Mn + 2) 

II 


X /I + 


4QL 


vMn+2 


AD„CT, 
s H p i aV 

2 Vi n Wi 

4QLBg(Mn + 2) 


(3 - n)[| 


Lf-'in 


(3 - n)(J 




(B17) 


Simplifying equation (B17) by virtue of equations (B5) and (B12a) gives 


3 AP 


s,l-2 




R *1 


4QL + 2f l LT lVA 


Q/A £ 


= 2 S P ! + P 2 ) A s C p¥l ‘ V T 2 “ V 


1 + 


t 2 - T r 

T, 


. Mn+1 


1 T 2 - T 1 
x /- v -A— 1 

w T, 


Pf-jLT-L 


A 


D H (Mn + 2) (T 2 - T ± ) 


(3 - n)2 
A 


^ + T 2 - Tq^ 2 


2f ;L LT 1 


D H (Mn + 2)(T 2 - T x ) 


(3 


- n aJ 


»i( 


(B18) 


Let t = T /T ; then equation (B18) becomes 
£ 


d AP 


s ,l-2 


df- 


Q/A c 


™ ( 2f L - 

= 2 R T 1 ) 4QL + _ 

8 P 1 + P 2 ( A s C p D H T l V* - 1> 


1 A T Mn+l 


[-(t - 1)] 


2f x L(3 - n)| t 


Mn+2 


2^(3 - n)f 


D (Mn + 2)(t - 1) D (Mn + 2)(t - l) 

n it 


(B19) 


Simplifying equation (B19) gives 


28 



d AP 


s ,1-2 


it) 


= 2 5 Tl 


4QL 


2f L ? 

+ i " 


Q/a £ 


S P 1 + P 2 A s C p D H T l 


D. 


H 


- 1 + 


(5 - n)r 


(5 - n) 


<M “ + ' (Mn + 2 )(t - Dt" 1 


(B20) 


Combining equations (B12a) and (B20) yields 


d AP 


s ,1-2 


it) 


= 2 - 


T 


Q/a £ 


g P + P A 
1 2 


. i 2f 1 Li Mn+1 

1 w ; x _ i + 1 


D. 


X 


- 1 + 


H 


5 - n / t 1 

Mn + 2 \t - 1 ~ , \ Mn+1 

\ (t - 1)t ) 


(B21) 


The minimum of the curves of AP against w/A occurs at a point 

S p _L — CL 

where 3 AP S ^ 1 . 2 /^^) = Oj hence, the sum of the terms inside the braces of 
equation (B2l) must equal zero for a minimum to exist. 


For most gases , the viscosity varies with the temperature to powers of 
less than one; that is , (a = where M < 1.0. 

For hydrogen gas in the range from 150° to 1500° R, M is approximately 

0.72. 


If the flow in the heated passage is completely turbulent (such that 
n = 0.2; it can be shown that 3 AP g (w/A) | ^ ^ is always positive. 


For laminar flow throughout the 

L/D h Of 430, a AP B; 1 . 2 /a(*/A) | Q/As 

negative for laminar flow at a value 
(i.e., the value of 


value of t 


crit 


is dependent on the ratio L/D ^ the 
ature, and the flow rate. 


heated passage, n = 1.0. Assuming a ratio 
in equation (B2l) changes from positive to 

of t of approximately 4.0. The exact 
t at which 3 AP n 9 /3 (w/a) I = 0) 

’ 'q/a s 

he at -flux value, the inlet fluid temper - 
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APPENDIX C 


SAMPLE CALCULATION 

The following sample calculation (run 30 ) shows the procedure used to 
obtain the calculated values in table II: 


Inlet static pressure, P , psia 9.87 

S ,1 

Inlet hydrogen gas temperature, T-^, °R 140 

Hydrogen gas flow rate, w, lb/sec 20.91X10" 6 

Outlet hydrogen gas temperature, T^.? °R 1004 

Inlet to exit static pressure drop, psi 0.289 


(a) Calculation of Reynolds number at test section inlet 
Tube diameter 


D h = 0.116 in. = 9.667X10" 3 ft 


Tube flow area 


A = (jt/4) (0.116) 2 = 0.010565 in. 2 = 7.34X10 -5 ft 2 
Viscosity of normal hydrogen gas (ref. 15) at test section inlet 

= 2.3X10 -6 lb/ (ft) (sec) 


Inlet Reynolds number 


R e inlet _ 


wD- 

Ap. n 


'E = ( 20 . 91X10" 6 ) ( 9 .667X10 -3 ) 
( 7 .34X10 -5 ) ( 2.3X10”®) 


Inlet Reynolds number 

■^ e inlet = 1200 

(b) Calculation of average heat flux 
Surface area of tube 


A 


= tcD^L = h Pl 116 ^ 50 ) = 0.1265 ft 2 
ET 144 


Rate of heat transfer to fluid 


/"T 2 

= w / C dl= w(H 2 - %_) 
*'T 1 P 


30 



From reference 15 


H 2 = f ( T 2’ P 2^ = 3425 Btu / 1 ' b 
H 1 = P R = 573 Btu / 1B 

Rate of heat transfer to fluid 

Q = 20 . 91X10“ 6 (3425-573) = 5.955X10" 2 Btu/sec 
Average heat flux 

— = 5 • 955x10 2 = 0.472 Btu/ (sec) (ft 2 ) 

A g 0.1265 

(c) Calculation of Reynolds number at test section outlet 
Viscosity of normal hydrogen gas (ref. 15) at test section outlet 

H 2 = 9.1X10 -6 lb/ (ft) (sec) 


M-l 

Re , . = Re. n , — 

outlet inlet j. 


Re , , = 1200 

outlet 


2.3X10 


-6 


9.1X10 


-6 


= 300 


(d) Calculation of fluid temperature ratio 


Fluid temperature ratio 

T = — = = 7.17 

T-l 140 


(e) Calculation of pressure-drop parameter 


AP = ( 9 . 87 ) (0 . 289 ) = 2 .85 (psia) (psi ) 

1 s.1-2 
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